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Abstract 

Many health outcomes are influenced by a person's body mass index, as well as by the trajectory of body mass index 
through a lifetime. Although previous research has established that body mass index related traits are influenced by 
genetics, the relationship between these traits and genetics has not been well characterized in people of South Asian 
ancestry. To begin to characterize this relationship, we analyzed the association between common genetic variation and five 
phenotypes related to body mass index in a population-based sample of 5,354 Bangladeshi adults. We discovered a 
significant association between SNV rs347313 (intron of NOSl AP) and change in body mass index in women over two years. 
In a linear mixed-model, the G allele was associated with an increase of 0.25 kg/m^ in body mass index over two years (p- 
value of 2.3-10"^). We also estimated the heritability of these phenotypes from our genotype data. We found significant 
estimates of heritability for all of the body mass index-related phenotypes. Our study evaluated the genetic determinants of 
body mass index related phenotypes for the first time in South Asians. The results suggest that these phenotypes are 
heritable and some of this heritability is driven by variation that differs from those previously reported. We also provide 
evidence that the genetic etiology of body mass index related traits may differ by ancestry, sex, and environment, and 
consequently that these factors should be considered when assessing the genetic determinants of the risk of body mass 
index-related disease. 
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introduction 

Body mass index (BMI) can predict subsequent healtli 
outcomes. Both over- and underweight individuals increase their 
subsequent risk of death and poor health outcomes when 
compared to their normal-weight peers [1-6]. Similarly, a change 
in BMI over time is also associated with increased mortality and 
morbidity [7-1 1]. BMI has increased in low-income countries over 
the last three decades [12], which has increased the percentage of 
these populations who are now exposed to the risks associated with 
increased BMI. Despite this increase in average BMI, much of the 
population in Bangladesh and South Asia remains underweight 
(including 40% of our population-based sample). Underweight 
individuals in Bangladesh also show increased mortality [4], and 
studies in other Asian cohorts have demonstrated that low BMI is 
associated with increased mortality from a wide variety of causes, 
including cardiovascular disease, cancer, and respiratory disease 
[13,14]. 



Although a person's BMI is affected by nutrition, genetics also 
exerts an influence on BMI-related phenotypes. This genetic 
component has been confirmed for point-in-time measurements of 
BMI related traits. Approximately 40 loci have been identified 
through genome-wide association studies (GWASs) [15-21], and 
heritability studies estimate that between 60%-80% of the 
phenotypic variance in BMI can be explained by genetic variation 
[22,23]. Although these conclusions come from studies that were 
predominantly composed of populations of European descent, 
GWASs of participants with African [24,25] and East Asian [26] 
ancestry suggest that associations between genetics and BMI- 
related traits can be found by studying a non-white population, 
such as our study sample. 

The Bangladeshi participants in this study can provide novel 
insights into the genetic mechanisms of BMI-related traits. The 
participants of this study have a genetic architecture unique to 
South Asia, and were also exposed to widespread undernourish- 
ment (one quarter of our participants had a BMI lower than 
17.6 kg/m ), both of which differentiates them from previously 
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studied populations. Their genetics and nutritional habits create 
the conditions by which some of the genetic determinants of BMI 
related traits may differ from those previously identified. This 
study also investigates a phenotype, change in BMI, which has not 
been widely examined with GWASs or heritability estimates. 

We conducted GWASs using genotyped and imputed data from 
a sample of Bangladeshi adults. We evaluated associations between 
genome-wide single nucleotide variants (SNVs) and the baseline 
phenotypes of BMI and height, as well as over- and underweight 
status. We also im'estigated whether there is evidence of a genetic 
component to the two-year change of BMI within the cohort. We 
repeated the GWAS analysis after we restricted the sample by sex 
and BMI classification at basefine, to investigate whether there was 
heterogeneity in the genetic determinants of BMI-related pheno- 
types. In addition, for each of the continuous phenotypes, we 
provided an estimate for the narrow sense (additive) heritability, as 
well as the phenotypic variability that is modeled by the additive 
effects of all measured and imputed SNVs. 

This study will help to increase the knowledge of the genetic 
determinants of BMI-related traits. We also introduce a BMI- 
related trait that has previously not been investigated as heritable: 
change in BMI, and suggest that genetic variation may drive this 
phenotype. We characterize the variants that contribute to BMI- 
related traits for the first time in South Asians, and suggest that the 
genetic variants that determine BMI related traits may differ than 
the variants identified in previously studied populations, and differ 
by gender. These results will help further elucidate the differences 
and similarities among populations in terms of the etiology of 
BMI. This biological insight may suggest how to control BMI- 
related morbidity, how to classify patients by risk of developing 
BMI-related diseases, and suggest populations that may be helped 
by interventions that may stabilize a person's BMI. 

Materials and Methods 

Study Sample 

The Health Effects of Arsenic Longitudinal Study (HEALS) 
[27] is a population-based cohort study established in 2000 in 
Araihazar, Bangladesh to primarily study the health effects of 
long-term arsenic exposure. The study enrolled a total of 20,033 
participants in two recruitment cycles. The Bangladesh Vitamin E 
and Selenium Trial (BEST) [28] was established in Araihazar and 
surrounding areas in 2006 to determine whether 7,000 arsenic- 
exposed participants who supplemented their diet with vitamin E 
and selenium would reduce their rate of non-melanoma skin 
cancer. Clinical and anthropomorphic measurements were ascer- 
tained biennially by trained clinicians and intervie\vc'rs, providing 
longitudinal information on the participants. Demographic and 
clinical characteristics of the study participants included in this 
analysis are shown in Table SI. 

Genotype Data 

A total of 5,499 participants from HEALS and BEST were 
genotyped using the lUumina Infinium: HumanCytoSNP- 1 2 v.2.1 
chip. Details of the study sample selection, quality control, and 
genotyping process are reported elsewhere [29]. The participants 
were genotyped in two batches. After quality control, 5,354 out of 
5,499 samples and 257,747 out of 299,140 SNVs were carried 
forward into imputation. 

MaCH software [30] pre-phased the genotyped data, and 
minimac software [31] imputed the data, using the HapMap3 
GIH reference panel. SNVs with a maximum likelihood r^ greater 
than 0.3 were used in analysis. A unix script assigned a hard call 
genotype probabilistically based on the minimac maximum 



likelihood probability. After imputation, 1,208,102 SNVs were 
used for association analyses. 

Assessment of Outcomes and Covariates 

At baseline and biennial follow-up, study physicians measured a 
participant's weight and height. Weight was measured three times 
over the course of the inter\dew with one of two Misaki scales that 
were manufactured in Japan (serial numbers: 67117, 58216). The 
scales were calibrated weekly. The average value of the three 
measurements was recorded. Height was measured by placing the 
scale on the participant's head, parallel to the floor, and measuring 
the length from the scale to the floor with a locally manufactured 
tape measure three times. The average value of the three 
measurements was recorded. BMI was calculated by dividing 
average weight in kilograms divided by the square of the average 
height in meters. Participants with a BMI less than 18.5 kg/m^ 
were classified as underweight, and participants with BMI greater 
than or equal to 23 kg/m^ were classified as overweight. A cutoff 
of 23 kg/ m^ was chosen for classification of overweight status (as 
opposed to 25 kg/ m^) based on evidence suggesting that people of 
Asian descent have higher adiposity and risk of obesity-related 
comorbidities at lower BMI than do people of European descent 
[32,33]. Gender and age were ascertained during the interview. 

Exclusions 

Among the 5,354 participants with GWAS data, baseline BMI 
was missing on 1 1 participants, baseline height was missing on 1 0 
participants, 103 participants were lost to follow up after baseline, 
and an additional 263 participants were missing a BMI 
measurement at their first follow up. Participants with missing 
information on a phenotype were excluded from analysis for that 
phenotype. 

Statistical Methods 

We analyzed five phenotypes: BMI, height, underweight and 
overweight (both were only compared to their normal weight 
peers), and change in BMI over two years. Previous studies suggest 
that the genetic underpinnings of BMI-related traits may differ by 
sex [35], so we also conducted subsample analyses on BMI and 
change in BMI after w(; stratified by sex. We also stratified by the 
participant's BMI at baseline (those who were underweight at 
baseline, those who were normal weight at basefine, and those who 
were overweight at baseline), and re-analyzed BMI and change in 
BMI for the participants in each category. 

Our interviews with the participants established that many 
participants were known to be distantly related, and our analysis of 
relatedness indicated that more than 60% of the participants were 
genetically related to at least one other participant as a third cousin 
or closer. Therefore, we used the software Genome-wide Efficient 
Mixed-Model Association (GEMMA) [34] to control for this 
between-subject genetic correlation. This software calculated a 
relatedness matrix based on the pairwise covariance between 
genotypes, and then estimated the effect of each SNV on the 
phenotype while controlling for the relatedness matrix with a 
Unear mixed model. We assumed an additive model of inheri- 
tance, and treated all phenotypes as continuous variables. We also 
controlled for the linear effects of baseline age, the square of 
baseline age, sex, and genotyping batch. The reported p-values are 
from a Wald test. We considered 5-10~" to be the cutoff for 
genome wide significance. R software [35] and LocusZoom [36] 
were used to plot results, and ANNOVAR [37] was used to 
annotate the context and nearest genes of the variants. 

We did not transform the continuous variables, as they are 
distributed roughly normally in this population, and the linear 



PLOS ONE I www.plosone.org 



2 



August 2014 I Volume 9 | Issue 8 | e105062 



GWASs and Heritability Estimates of BMI-Related Traits in Bangladesh 



mixed model technique used by GEMMA is relatively robust to 
deviations from normality [38]. 

For the continuous traits, we estimated three aspects of 
heritability (described below) using Genome-wide Complex Trait 
Analysis software (GCTA) [39]. Similar to GEMMA, this software 
also estimates a relatedness matrix based on the pairwise genetic 
covariance. We first estimated hg^, the amount of variance in the 
trait that was explained by the interrogated SNVs. For hg^, GCTA 
fits the given phenotype with a linear mixed model, while it uses 
the estimated relatedness matrix as the variance term, to estimate 
the variance explained by all SNVs. We repeated the hg^ analysis 
after we restricted our sample to only one participant of any pair 
where the estimated kinship coefficient was larger than 0.025 
(approximately 2"'^-3'^'' cousins), to reduce the possibility that the 
estimates would be inflated because of shared environment. We 
also estimated the full narrow sense heritability, h^, using the 
method described by Zaitlen et al [40] . This procedure is similar to 
the one described above, but replaces the fuU relatedness matrix 
with a modified one that assumes zero relatedness between 
participants whose estimated relatedness is less than 0.05. This 
modified relatedness matrix better approximates the identity by 
descent matrix. 

Ethics Statement 

The study procedures and consent procedures were approved 
by the Columbia University Institutional Review Board and the 
Ethical Committee of the Bangladesh Medical Research Council. 
Since many of the people in rural Bangladesh are unable to read, 
each individual who agreed to participate provided verbal consent, 
which was recorded by the field staff physicians in the interview 
form in the presence of a witness. The study team explained details 
of the study procedures and the benefits and risks of the study in 
local language. Participants were advised that they could consent 
with or without donating blood or urine. The study team also 
explained to the participants that they could withdraw from the 
study at any stage, even if they had already provided consent, and 
explained the procedure for withdrawal. 

Results 

Genome-Wide Association Analyses: Baseline BMI and 
Height, and Baseline Over- or Under-weight Status 

We carried out GWASs on four BMI-related phenotypes that 
were measured at the baseline visit: BMI, height, overweight status 
and underweight status. No SNVs were associated with any of the 
baseline traits with a p-value that reached the genome-wide 
threshold. Table 1 lists SNVs associated with all phenotypes for 
which we estimated a p-value of less than 10 QuantUe-quantile 
plots suggest that the p-values for baseline BMI are enriched for 
small p-values, compared to the null hypothesis of no genetic 
association (Figure SI a). 

Ten previously published GWASs investigated cross-sectional 
BMI and reported a genome-wide significant association for a 
variant within or around the fat mass and obesity-associated 
protein (FTO) gene [13-15,21,22,36-40]. Visual inspection of our 
data does not show convincing evidence for a strong signal in the 
FTO region (Figure 1). 

Our BMI GWAS interrogated forty-seven SNVs that the 
National Human Genome Research Institute (NHGRI) GWAS 
catalogue lists as associated with BMI with a p-value smaller than 
10 ^. Thirty-four of these overlapping SNVs (72%) have an 
estimated effect size in our analysis that is directionally consistent 
with previous reports. While we estimated a nominally significant 
(p<0.05) p-value for only 9 of these overlapping SNVs, the SNVs 



whose estimated effects are in the same direction as previous 
literature have significance tests that are enriched for small p- 
values (Figure 2A). There is no such enrichment in the quantile- 
quantile plot of the SNVs where our reported effect direction is 
different than the reported literature (Figure 2B). 

We next stratified the sample by sex and re-analyzed the 
GWAS. No SNVs reached the genome-wide threshold (Table 2). 
QuantUe-quantile plots of the stratified analysis suggest an inflation 
of small p-values for baseline BMI for females (Figure S2b), and 
baseline overweight status for males (Figure S2g). 

Lastly, we examined whether the association between common 
variation and BMI differed by baseline BMI category. In these 
restricted GWAS, the analysis of the 869 participants who were 
overweight at baseline identified three SNVs that were genome- 
wide significant, and eight more that were suggestive of genome 
wide significance (p<10 ^) (Table S2, Figure S2c). However, after 
closer examination, the participant in our study with the largest 
measured BMI carried the risk allele at each of these loci. Even 
though the BMI measurement for this participant seemed to be 
accurate, we were concerned that this one observation was driving 
the observed association (this participant's BMI was 51 kg/m^, 
and the participant with the second largest BMI measured 37 kg/ 
m^). When we repeated the analysis excluding this one participant, 
no SNV rose to the level of genome-wide significance (Table S2, 
Figure S3d). 

Genome-Wide Association Analysis: Change in BIVll 

On average, our study participants increased their BMI 
0.204 kg/m^ (confidence interval: 0.17-0.24 kg/m^) between their 
initial interview and their first follow up two years later. At 
baseline, increasing age was associated with a small decrease in 
BMI (Pearson correlation coefficient: -0.063, p-value: <0.0001), 
which suggested that the within-person increase was not a result of 
age-related weight gain. 

We assessed with a genome-wide association study whether 
genetic variation affected a person's propensity to change their 
BMI over time. When examining the entire sample, no SNV was 
associated with this change in BMI phenotype with a p-value that 
reached the level of genome wide significance, although three of 
the SNVs were suggestive of association (Table 1). 

As with the cross sectional BMI-related phenotypes, we 
investigated whether the SNVs that drove the change in BMI 
differed by sex, and found evidence that they do. The genotyped 
SNV rs347313 on chromosome 1 was significantly associated with 
two-year change in BMI at the genome-wide level (p = 2.3- lO"'') 
(Figure 3). This G/A SNV is in the intron of the gene nitric oxide 
synthase 1 adaptor protein (NOSIAP) (Figure 4). The risk allele, 
G, was associated with a per-aJIele increase of0.25kg/m2inBMI 
over the two-year period. This SNV was not associated with the 
cross-sectional measurements of baseline BMI, baseline weight, 
baseline height, or being over- or underweight at baseline in our 
analyses, and was also not associated with change in BMI for 
males (p>0.05 for all associations). 

As with the cross-sectional BMI phenotypes, we also investigat- 
ed whether the SNVs that are associated with a change in BMI 
differed by baseline weight, and found no evidence that they did 
(Table S2, Figure S3). 

Heritability 

We estimated heritability for each of the three continuous traits 
in three ways: we used two methods to estimate hg^, which is the 
upper limit of the amount of phenotypic variance we could have 
expected to explain with our GWAS; and we estimated h^, the fuU 
narrow sense heritability. 
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Figure 1 . P-values for Baseline BWII in tlie Vicinity of tlie FTO Gene, p-values in the vicinity of the FTO gene indicate no obvious association 

with body mass index near the gene. 

doi:10.1371/journal.pone.0105062.g001 
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Figure 2. QQ Plots for SNVs Previously Associated with Baseline BIVll. QQ plots for the 47 SNVs previously reported as associated with body 
mass index that we interrogated (A) SNVs that have a directionally consistent effect with previous literature (n = 35) show an inflation in the p-values 
calculated for this study, whereas the QQ plot of our p-values for SNVs that are directionally inconsistent (B) (n = 13) are consistent with the null 
hypothesis. 

doi:1 0.1 371/journal.pone.01 05062.g002 



PLCS ONE I www.plosone.org 



5 



August 2014 | Volume 9 | Issue 8 | el 05062 



GWASs and Heritability Estimates of BMI-Related Traits in Bangladesh 



> 

Q. 



< 



< 



> 

z 



Q_ ^ ^ 

I— M- M- 



^ ^ rf 

.t: rsi m .-hi 



ro ro ro ro 



CD Q) C7) c 
S U S _^ 



o ro 
o LT) Lo 

Cr\ fN 0^ ro 



vo o 



0^ LO (Ji LO LO 

L/1 r~v Lo r^_^ fN 
r-j fN r-j rs rN 



t— LO 1— 



ro ro ro ro 



£ ^ 



— E 



5 g 



£ E £ E 



oi O) ai ai ■- -- 



g g s s 



CT) cn CTi 01 



O O O O U U 



5 g 



£ i5 



a; £ 



Al 
1 

11 

(D 

i o 

A "J 

£ -o 

o c 

^ . . 

u E 

5 I 

£ o 



t| 

ii 

o 



OJ 



Q. £1 



■ft, < 



"D O 
OJ c 

OJ rs 

I ■£ S 

5i 5 

" ^ o 

^ Q. 

^ c 75 



OJ 



■Sg^ 

i_ OJ ro 

^ ■§? 

< !S o 



PLOS ONE I www.plosone.org 



6 



August 2014 | Volume 9 | Issue 8 | el 05062 



GWASs and Heritability Estimates of BMI-Related Traits in Bangladesh 




Figure 3. Manhattan and QQ Plots of Change in BMI Over 2 Years. A signal at chromosome 1 , and inflation in the tails of the QQ plot indicate 
an association between SNV rs347313 and change in BIVll over 2 years in females. This SNV is in the intron of N0S1AP. The suggestive threshold is 
drawn at 10"^ and the genome wide significance threshold is drawn at 5-10"^. 
doi:1 0.1 371 /journal.pone.01 05062.g003 
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doi:1 0.1 371/journal.pone.01 05062.g004 



PLOS ONE I www.plosone.org 



7 



August 2014 | Volume 9 | Issue 8 | el 05062 



GWASs and Heritability Estimates of BMI-Related Traits in Bangladesh 



We found evidence that the SNVs interrogated by our 
genotyping and imputation were associated witli tlie phenotypic 
variation in the continuous BMI-related traits of this study 
population (Table 3, estimates). We repeated the hg"analysis, 
including one participant of any pair where the estimated kinship 
coefficient was larger than 0.025 to remove inflation in the 
hg^estimates that could be a result of phenotypic similarity due to 
shared environment. While this analysis reduced our sample size 
by half, the estimates were consistent with our full-sample hg , 
albeit with wider standard errors. 

We estimated the full narrow sense heritability using the method 
described by Zaitien et al (Table 3, h^ estimates). The heritability 
estimates from that method are strikingly high, suggesting that the 
narrow sense heritability is dominated by shared environment, 
epistatic interactions, or dominance effects, and that this method 
might not be an appropriate estimate of the narrow sense 
heritability in our population. 

As with the GWAS analyses, we repeated the heritability 
analysis, separating the participants by sex (Table 3). For all traits, 
men tended to have more phenotypic variance explained by the 
SNVs, although the estimates were imprecise. We also repeated 
the analysis, separating the participants by their baseline BMI 
categorization (Table S3). 



3 ■?! 



Discussion 



In this analysis, we observed suggestive associations between 
several SNVs and BMI-related phenotypes measured at the 
baseline of our study. This begins to characterize the genetic 
determinants of BMI-related phenotypes for the first time in a 
South Asian population, although most associations did not reach 
genome-wide significance, and multiple phenotypes were tested. 
For some of the traits, the quantUe-quantile plots are inflated, and 
the heritability estimates for the continuous phenotypes are 
significant, which suggests that a larger sample size may be able 
to detect SNVs with modest effect. 

The results of our stratified analysis also suggest that the genetic 
drivers of BMI-related traits may vary by sex, and BMI 
categorization. For each of the sub-analyses, the SNVs that our 
stratified analyses identified as suggestive (p<10~^) did not overlap 
with the SNVs identified as suggestive in the pooled analysis. Since 
in the Bangladeshi population, both over- and underweight 
individuals increase their risk of morbidity from both communi- 
cable and non-communicable diseases, these results may help to 
develop risk scores that identify segments of a population that may 
be most susceptible to illness, allowing for more targeted 
intervention of anti-obesity and anti-malnourishment therapies. 

These findings will be informative about which genetic 
determinants of BMI are consistent across populations, and which 
differ. It is compelling that our investigations show that many of 
the loci previously indicated as being associated with BMI [15,17- 
20,41,42] contain no SNVs that are strongly associated with BMI 
in our sample. While the participants of these studies were 
primarily of recent European descent, one study in an African 
American population [21] and one in a Japanese population [22] 
found an association signal from a variant in the vicinity of FTO, 
which replicated a signal seen in populations of European ancestry 
[15,18,19,41,42]. This suggests that genetic variation at this locus 
affects BMI in people of difiFering backgrounds. These cUfferences 
between the previous research and our study may be a result of the 
mechanisms by which genetics influences BMI. There is growing 
evidence that many of the variants identified in previous BMI 
GWASs affect BMI by influencing food choices or appetite 
[19,43,44]. However, in the rural Bangladeshi population. 
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decisions about food may be driven less by preference, and more 
by availability and alfordabiUty, limiting the ability of any SNV 
that modifies behavior to affect BMI. 

This analysis also reports the first genetic variation significandy 
associated with change in BMI that was discovered by interrogat- 
ing the whole genome. The change in BMI for women with the 
risk allele is modest, 0.25 kg/m^ over two years, but weight loss as 
low as 0.2 kg/year, and gains as little as 0.6 kg/year have been 
found to be associated with all-cause mortality in a different 
population [45]. To our knowledge only one previous GWAS 
investigated adult change in BMI; it reported a null result [15]. 
The results presented here represent a much denser array of 
SNVs, and also have benefited from improvements in analysis 
since that study was published, such as control for population 
substructure, and imputation of non-genot^ped SNVs. As this is 
one of the first studies to examine change in BMI in a genome- 
wide way, these findings need to be replicated. If the signal is 
replicated, further research needs to examine whether this same 
variant is also responsible for a genetic contribution to change in 
BMI in women in other populations, or whether this variant affects 
females" BMI trajcctorv' only in the context of the Bangladeshi 
environmental or genetic background. 

The SNV we identified, rs347313, has never been found to be 
associated with a complex trait phenotype, according to the 
National Human Genome Research Institute GWAS catalogue. 
This SNV is situated in the intron of the gene NOSIAP (also 
known as CAPON) [46], which encodes a cytosolic protein that 
binds to neuronal nitric oxide synthase, a signaling molecule. 
NOSIAP has been associated with cardiac phenotypes [47-51], 
and marginally associated with childhood hip circumference in a 
Hispanic population (p = 8.6-10 ^ [52]). 

Our weight measurement protocol likely captured some intra- 
individual variability in weight, since participants were interviewed 
at different times of the day, and during difierent seasons of the 
year. However, this variation in the timing of the weight 
measurement is likely to be smaller than the inter-individual 
variation in weight, and also not associated with genetics. 
Therefore, we expect that the noise introduced by this weight 
measuring protocol would not be systematic, and would not result 
in confounding between genetics and body mass index. A previous 
study of BMI and mortality on a subset of the same study 
population was associated with other clinical endpoints in the 
expected manner [4] . Therefore we have confidence that the BMI 
we measured accurately captures the BMI of this population. 

In future research, we hope to extend this investigation to 
examine whether genetics interacts with nutrition in such a way 
that different variants affect BMI in the presence of different 
nutrition choices. In the near future, we wiU have available data on 
the food intake of these participants, and we intend to investigate 
the possibility that nutrition mediates the influence of genetic 
variation on BMI, and the possibility of nutrition interacts with 
genetic variation to influence BMI. 
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